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Summary 


Bismuth  telluride  and.  its  alloys  are  the  most  important  thermo¬ 
electric  materials  for  applications  near  room  temperature.  In  the  past, 
it  has  always  been  made  n-type  by  adding  a  foreign  impurity.  In  the 
first  section  of  this  report,  some  measurements  are  reported  on  material 
made  n-type  by  adding  excess  tellurium.  At  very  low  temperatures,  the 
electron  mobility  is  about,  ten  times  larger  than  has  been  previously 
reported.  In  the  temperature  range  where  the  material  is  a  useful  thermo¬ 
electric  material,  there  is  little  change  in  the  mobility.  However,  the 
larger  mobility  at  helium  temperatures  may  allow  one  to  do  experiments 
such  as  cyclotron  resonance  which  will  give  a  better  measurement  of  the 
electronic  band  structure  than  has  been  possible  previously. 

The  second  paper  gives  the  results  of  an  investigation  of 
the  Seebeck  coefficient,  thermal  conductivity,  and  electrical  conductivity 
of  bismuth  single  crystals.  Bismuth-antimony  alloys  have  been  used  as 
thermoelectric  device  elements  for  many  years,  and  they  have  the  highest 
figure  of  merit  of  any  material  at  temperatures  near  100°K.  In  this 
work,  some  of  the  properties  of  the  band  structure  and  the  lattice 
thermal  conductivity  have  been  determined.  The  optimum  figure  of 


merit  of  the  material  has  been  calculated. 
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HIGH  MOBILITY  N-TYPE  BISMUTH  TELLURITE 

Roland  W.  Ure,  Jr. 

Westinghouse  Research  Laboratories 
Pittsburgh  35,  Pa.,  U.S.A. 

ABSTRACT 

Samples  of  n-type  bismuth  telluride  with  carrier  concentrations 
between  3  x  lO1^  and  6  x  lO1^  cm'^  have  been  made  by  doping  with  excess 
tellurium.  The  largest  mobility  at  4.2°K  iB  1.8  x  10^  cm2/Vsec,  approxi¬ 
mately  ten  times  larger  than  has  been  observed  in  n-type  iodine  doped 
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bismuth  telluride.  At  a  carrier  concentration  of  about  5  ,x  10  cm0, 
there  is  a  sharp  change  In  slope  of  the  curve  of  4.2°K  mobility  vs 
carrier  concentration.  Th.se  results  are  interpreted  on  the  basis  of 
a  two-band  model.  Seebeck  coefficient  measurements  tend  to  confirm  this 
model. 


I.  INTRODUCTION 

The  compound  BigTe^  1b  an  anisotropic  semiconductor  with  a 
rhombohedral  structure  and  a  very  pronounced  cleavage.  When  grown  from 
a  stoichiometric  melt,  it  is  p-type  with  a  carrier  concentration  of  the 
order  of  lO1^  holes /cm^.  The  reason  for  this  high  carrier  concentration 
is  that  the  maximum  melting  composition  is  on  the  bismuth  rich  side  of  the 
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stoichiometric  composition  (offergeld  and  Van  Cakeriberghe,  1959)  and  hence 
a  crystal  grown  from  a  stoichiometric  melt  contains  excess  Bi.  This  excess 
Bi  aats  as  an  acceptor  in  the  crystal.  Measurements  of  a  number  of  different 
properties  of  both  n  and  p-type  saterial  have  been  made  (for  review  of  the 
properties  of  Bi^Te^,  see  Ure;  1961)  but  the  carrier  concentration  is 
usually  changed  by  adding  a  foreign  impurity  such  as  iodine.  We  have 
reported  earlier  some  data  on  n  and  p-type  material  in  which  the  carrier 
concentration  was  altered  by  changing  the  Bi:Te  ratio  in  the  cryBtal 
(Satterthvaite mid  Ure,  1957).  Sladek  (1959)  has  measured  one  of  our 
samples  and  found  a  4.2°K  mobility  approximately  ten  times  larger  than  has 
been  reported  elsewhere. 

In  this  paper,  we  report  measurements  of  Hall  coefficient  and 
mobility  over  a  wider  range  of  temperature  (U.2°K  to  250°K)  and  carrier 

concentration  (3  x  lO1^  to  6  x  lO1^)  tor  n-type  samples.  We  find  mobilities 
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as  high  as  180,000  cm  /Vsec  at  4.2JK.  At  room  temperature,  the  mobility  is 
about  the  same  as  for  iodine  doped  samples. 

II.  EXPERIMENTAL 

All  samples  measured  were  single  crystals  and  were  grown  from 
melts  containing  from  63  to  67  atomic  percent  tellurium.  Both  Bridgman 
and  zone  leveling  techniques  were  used.  Considerable  difficulty  was 
experienced  in  obtaining  reasonably  homogeneous  samples,  particularly 
at  carrier  concentrations  below  101®  cm"^.  In  samples  grown  by  the 
Bridgman  procedure,  there  is  a  regular  increase  in  the  carrier  concen¬ 
tration  from  die  bottom  to  the  top  of  the  crystal  because  of  the  increase 

in  Te  content  of  the  liquid  as  the  crystal  freezes.  However,  some  uniform 
samples  were  obtained  bv  cutting  samples  perpendicular  to  the  ecrowth  direction 


ware  'tut  with  ft  apart;  cutter  (Chandrasekhar,  1961} 

Cole,  Bucklov,  and  Gr?.  pyon ,  1963,).  The  edgec  of  the  samples  vhich  were 
damaged  during  the  noting  opejratioi:  -rere  cleaved  off  and  the  samples 
were  etched  in  Ini  EOItjEE-TO.  tc  class  the  surfaces  and  to  detect  cracks 
in  the  samples.  Cmc'A.~:?V',t.  aamj-l eo  were  usually  obtained  by  this  pro¬ 
cedure.  thwaple  irwensicv,*  wsr-.t/J.  iron  about  1  x  1  x  1C  ma  to  2  x  3  x  20  mm. 

The  er-tfe  of  the  cemples  were  nickel  plated  end  platinum  current 
leads  were  soldered  to  the  plating.  1 let i nun- rhodium  pressure  contacts 
were  used  for  potential  p robin .  Four  probes  ware  placed  on  the  samples  — 
two  on  each  eide  :\o  that  two  resistivity  and,  two  Hall  measurements 
could  be  made.  On  ch«  bettor  samples  (indicated  by  filJLed  points  on 
the  curves)  the  two  Hall  measurements  differed  by  less  than  1 %  and  the 
two  resistivity  waafiuvinarnbo  differed  by  less  than  10J&.  All  measurements 
were  nuide  by  Htavidurd  do  pot  st*  biometric  methods.-  lr,  measuring  resistivities, 
the  readings  were  first  made  vil'-r  the  current  flowing  in  one  direction. 

Then  a  reading  wan  taken  immediately  after  reversing  the  current  before 
the  temperature  distribution  in  the  r-wnple  changed.  In  this  way,  thermal 
emf ’s  resulting  from  the  temperature  Inhoaogeneity  produced  by  the  Peltier 
effect  were  eliminated.  Hall  effect  measurements  vere  made  in  a  field  of 
6  kG.  Readings  were  taken  with  both  currant  and  field  in  the  forward  and 
reverse  directions,  and  the  appropriate  average  of  four  readings  was  used. 

All  measurements  were  wade  with  the  magnetic  field  perpendicular  to  the 
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cleavage  planes  and  the  current  in  the  cleavage  plane.  Temperatures  were 
measured  with  a  carton  resistor  or  with  a  copper-constantan  thermocouple 
both  calibrated  against  a  standard  platinum  resistor. 

The  measured  Hall  coefficients  are  shown  in  Jig.  1.  Carrier 
concentrations,  n,  were  calculated  from  the  equation 

n  *  Br/Re  (l) 

using  the  average  of  the  two  Hall  measurements  R  at  77°K.  Here  B  is  a 
factor  depending  on  the  band  structure  and  is  unity  for  spherical  bands. 
From  magnetoreaistance  measurements  on  iodine  doped  samples,  Drabble, 
Groves,  and  Wolfe  (1958)  showed  that  the  BigTe^  has  a  multivalley  con¬ 
duction  band  with  highly  oblate  constant  energy  surfaces.  They  found 
B  =  0.326.  However,  Goldsmid  (1961)  and  Walker  (i960)  have  found  that 
this  factor  is  a  function  of  the  carrier  concentration  in  the  specimen. 
Because  of  the  significant  differences  between  the  properties  of  iodine 
doped  and  tellurium  doped  samples  reported  here,  the  factor  may  be 
different  in  our  samples.  Because  of  these  uncertainties,  we  have 
not  corrected  our  data,  but  have  taken  B  «  1.  The  factor  r  depends  on 
the  scattering  mechanism  and  the  degree  of  degeneracy.  We  have  taken 
r  «  1  since  the  scattering  mechanisms  have  not  been  established.  Ball 
mobilities  p  were  calculated  from  u  *  R^/p.  Values  of  the  Hall  mobility 
are  shown  in  Fig.  2  and  3  and  the  resistivity  p  in  Fig.  U.  It  must  be 
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remembered  that  the  "carrier  concentration"  quoted  may  be  3  times 
larger  than  the  true  carrier  concentration  and  that  the  Hall  mobility 
may  be  one-third  of  the  conductivity  mobility. 

III.  DISCUSSION 

A.  Hall  Coefficient 

The  Hall  coefficient,  shown  in  Fig.  1,  ehowa  some  variation 
with  temperature  in  the  extrinsic  region.  With  increasing  temperature, 
the  Hall  coefficient  first  falls  slightly,  goes  through  a  shallow  minimum 
near  the  temperature  of  liquid  nitrogen,  and  then  goes  through  a  slight 
maximum  before  it  falls  rapidly  as  the  samples  become  intrinsic.  Yates 
(1959)  observed  a  similar  behavior  in  iodine  doped  bismuth  telluride. 

In  a  semiconductor  which  has  significant  carrier  densities  in  each  of 
two  different  conduction  bands  with  extrema  at  different  energies  and 
with  different  mobilities,  the  Hall  coefficient  vs  temperature  shows  a 
maximum.  Thus  the  observed  maximum  is  evidence  for  the  two  band  picture 
to  be  discussed  later.  The  fact  that  sample  D-13  with  the  lowest  carrier 
concentration  does  not  show  this  maximum  may  be  evidence  that  the  energy 
separation  between  the  two  bands  is  large  enough  that  there  are  very  few 
carriers  in  the  upper  band  in  this  sample.  However,  this  sample  was 
fairly  inhomogeneous. 
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B.  Hall  Mobility 

The  largest  Hall  mobility  observed  is  1.8  x  10^  cm^/Vsec  in 
sample  B-2A  with  a  "carrier  concentration"  of  9  x  10^  cnA  This  sample 
was  reasonably  uniform.  The  two  Hall  coefficients  differed  by  7  percent 
and  the  two  resistivity  measurements  by  one  percent.  This  Hall  mobility 
is  ten  times  larger  than  the  largest  mobility  observed  by  Yates  (1959) 
on  iodine  doped  samples. 

In  the  region  between  77°K  and  200°K,  the  mobility  varies  as 
T~a  where  a  is  a  constant.  The  values  for  a  are  given  in  Table  1. 
Drabble,  Groves,  and  Wolfe  (1958),  Goldsmid  (1958),  and  Mansfield  and 
Williams  (1958)  have  all  found  values  very  close  to  T  for  iodine 

doped  material. 

The  most  obvious  explanation  for  these  differences  is  that 
iodine  doped  samples  still  have  an  appreciable  density  of  excess  bismuth 
acceptors  and  hence  are  heavily  compensated.  Thus  in  an  iodine  doped 
sample  with  a  carrier  concentration  of  10  cm  ,  there  may  be  of  the 
order  of  10^  ionized  bismuth  acceptors  plus  1.1  x  10^  ionized  iodine 
donors.  Thus  the  density  of  ionized  impurity  centers  is  of  the  order 
of  20  times  as  large  as  the  carrier  concentration.  This  ratio  of  the 
density  of  ionized  impurities  to  the  carrier  concentration  will  decrease 
as  the  carrier  concentration  is  increased  by  adding  additional  iodine. 

In  the  tellurium  doped  samples  this  ratio  is  one.  Thus  we  expect  that 
the  difference  between  the  mobility  in  iodine  and  tellurium  doped  samples 
will  be  largest  for  low  carrier  concentrations  and  will  decrease  as  the 
carrier  concentration  increases.  As  shown  on  Fig.  3  this  is  indeed  the 


case. 
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Since  ionized  impurity  scattering  becomes  less  Important  with 
Increasing  temperature,  the  difference  In  the  temperature  variation  of 
the  mobility  in  iodine  and  tellurium  doped  samples  can  be  explained  by 
assuming  that  there  is  still  a  small  amount  of  ionized  impurity  scattering 
at  100°K  in  the  compensated  iodine  doped  samples .  This  idea  is  supported 
by  the  fact  that  the  mobility  at  77°K  for  iodine  doped  samples  is  25  to 
35  percent  smaller  than  in  our  tellurium  doped  samples.  [Data  on 
and  resistivity  at  77°K  in  iodine  doped  samples  taken  from  Drabble, 

Grover,  and  Wolfe  (1958)  and  Bovley,  Delves,  and  Goldsmid  (1958) ]• 

However,  Bowley,  Delves,  and  Goldsmid  (1958)  have  determined  the  value 
of  s  in  the  expression  for  the  relaxation  time  t  =  tqE8  from  measure¬ 
ments  of  the  magneto-8eebeck  effect,  and  they  find  a  value  close  to  the 
theoretical  value  for  acoustic  mode  lattice  scattering.  The  difference 
between  their  measured  value  and  the  theoretical  value  indicates  that, 
although  the  amount  of  ionized  impurity  scattering  is  small,  the  amount 
of  ionized  impurity  scattering  increases  with  decreasing  carrier  con¬ 
centration.  Austin's  (i960)  free  carrier  absorption  data  indicate 
some  ionized  impurity  scattering  at  77°K  in  iodine  doped  samples. 

C.  Magnitude  of  Mobility  at  U.2°K 
In  the  degenerate  case,  the  mobility  due  to  ionized  impurity 
scattering  (Brooks,  1955;  Dingle,  1955}  and  Mansfield,  1956)  is 


8 


Hj  =  3€^3n/l6n2amV3l,if(l))>  (2) 

where 

f(b)  =  -£n(l+b)-'b(l+h)"1, 

b  =  (3n/n)l/3(eoh2/4ode2), 

cQ  is  the  static  dielectric  constant,  n  is  the  carrier  concentration, 
m^  is  the  mobility  effective  mass,  md  is  the  density-of-states  effective 
mass,  and  is  the  concentration  of  scattering  impurities.  Brooks  (1955) 
has  discussed  the  relation  between  the  mobility  mass  and  the  mass  tensor 
for  the  individual  valleys.  For  the  purposes  of  our  estimates  we  will 
take  mm  =  md>  Austin  (1958)  has  measured  an  optical  dielectric  constant 
of  85.  Since  the  static  dielectric  constant  will  be  larger  than  this, 
we  will  assume  a  value  of  100.  In  order  to  have  a  mobility  of  1.8  x  10^ 
cm  /Vsec  with  10  carriers  an  effective  mass  of  0.055  ®Q  ®uat  be  used. 

In  Table  II  we  have  shown  values  of  the  effective  mass  derived  from 
Walker's  data  (i960)  on  iodine  doped  material  and  our  measurements  on 
one  sample  shown  in  Fig.  5.  Here  we  have  used  the  standard  theory  for 
the  Seebeck  coefficient  assuming  that  the  relaxation  time  is  of  the  form 
t  =  tqEb  .  Values  are  given  for  s  =  l/2,  the  value  for  acoustic  mode 
lattice  scattering,  and  s  =  3/2,  the  simple  theory  for  ionized  impurity 
scattering.  The  values  for  other  scattering  mechanisms  lie  between  these 
two.  Since  the  predominant  scattering  mechanism  in  our  sample  at  100°K 
is  probably  acoustic  mode  lattice  with  small  amounts  of  optic  mode 
scattering,  a  reasonable  estimate  of  the  density-of-states  mBBB  is  0.2  mQ. 
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With  this  mass,  a  dielectric  constant  of  360  must  be  assumed  in  order 
to  get  agreement  between  Eq.  2  and  our  experimental  values.  In  calculating 
the  effective  masses  given  in  Table  II,  the  constant  B  in  Eq.  1  was  taken 
as  1.  If  the  value  of  B  found  by  Drabble,  Groves,  and  Wolfe  (1958)  were 
used,  the  effective  masses  shown  in  Table  II  would  be  multiplied  V*  0.47* 
Howe-ir,  to  be  consistent,  we  have  to  fit  a  mobility  of  5*9  *  10  instead 
of  1.8  x  10^,  and  the  dielectric  constant  required  is  still  300.  Thus 
there  seems  to  be  three  possibilities:  (l)  Bismuth  telluride  has  a  static 
dielectric  constant  which  is  much  larger  than  its  optical  dielectric  constant, 

(2)  the  present  theory  of  ionized  impurity  mobility  is  not  accurate,  or 

(3)  our  estimates  of  the  effective  mass  are  too  large.  A  similar  situation 
has  been  observed  in  lead  telluride  (Allgaier  and  Scanlon,  1958)*  From 
measurements  of  the  polaron  minimum  in  tunneling  experiments  on  lead 
telluride.  Hall  (i960)  has  concluded  that  the  static  dielectric  constant 

is  not  much  larger  than  the  optical  dielectric  constant  in  this  material. 

A  number  of  objections  to  the  present  theory  of  ionized  impurity 
scattering  have  been  expressed  (Herring,  i960).  However,  in  our  case 
the  Born  approximation  should  be  valid  since  the  dielectric  constant 
is  100  or  larger. 


D.  Variation  of  4.2  K  Mobility  with  Carrier  Concentration 

As  shown  on  Fig.  3>  the  mobility  at  4.2°K  varies  with  carrier 

“O  28  10 

concentration  as  n"w*  for  "carrier  concentrations"  less  than  5  x  10 
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and  as  n  for  higher  "carrier  concentrations."  The  slope  at  the  lower 
carrier  concentration  agrees  with  the  theory  of  Eq.  2  within  the  experi¬ 
mental  scatter  of  the  data. 

A  rapid  decrease  in  motility  with  carrier  concentration,  similar 
l8 

to  our  data  for  n  >  5  x  10  ,  has  been  found  in  FbTe  (Kanai,  Nii,  and 

Watanabe,  1961).  It  has  been  explained  by  assuming  that  the  scattering 
cross  section  remains  independent  of  carrier  concentration.  The  theory 
of  ionized  impurity  scattering  shows  a  rapid  dependence  of  the  relaxation 
time  on  the  carrier  velocity  and  thus  a  rapid  variation  of  scattering 
cross  section  with  carrier  concentration  in  degenerate  materials. 
Vinogradova,  et  al.  (1959)  have  found  that  the  scattering  cross  section 
for  ionized  impurities  in  a  80jfc  BigTe^-20]k  BigSe^  alloy  at  temperatures 
of  100  to  300°K  is  independent  of  electron  energy.  They  explain  this  by 
assuming  that  the  large  dielectric  constant  screens  the  Coulomb  field 
sufficiently  so  that  the  major  scattering  is  from  the  core  of  the  impurity 
rather  than  from  the  Coulomb  field  of  the  Impurity.  This  cannot  be  the 
explanation  for  our  results  at  4.2°K  since  the  theoretical  scattering 
by  the  screened  Coulomb  field  alone  gives  a  mobility  smaller  than  that 
observed  experimentally. 

We  are  tempted  to  explain  the  change  in  slope  of  the  curve  of 

Fig.  3  as  a  two  band  effect.  For  carrier  concentrations  less  than  about 
l8 

5  x  10  cm  ,  all  of  the  electrons  are  in  one  band.  At  higher  concen¬ 
trations,  the  electrons  begin  to  fill  a  second  band  having  a  higher 
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effective  mass  and  lover  mobility.  As  donors  are  added  and  the  total 
carrier  concentration  increases,  the  density  of  scattering  centers 
becomes  larger  than  the  carrier  concentration  in  the  lower  band  n^, 
and  the  nobility  of  the  elec  -one  in  the  lower  band  decreases.  Also 
since  the  upper  band  is  assumed  to  haw  a  larger  mass,  the  electrons  in 
the  upper  band  will  have  a  lower  mobility  than  those  in  the  lower  band. 

Thus  -she  average  mobility  will  decrease  with  increasing  carrier  concentration. 
Preliminary  calculations  taking  into  account  the  proper  relation  between  the 
Hall  coefficient  and  the  carrier  concentration  in  each  band  seem  to  show 
that  the  rate  of  reduction  in  mobility  should  be  smaller  than  observed 
experimentally.  However,  additional  calculations  are  in  progress. 

Yates  (1959)  and  Walker  (i960)  have  suggested  that  a  two  band  model 
would  explain  some  of  the  features  of  their  work.  The  fact  that  the 
effective  masses  shown  in  Table  II  for  more  heavily  doped  material  are 
appreciably  larger  than  the  masses  calculated  for  our  lightly  doped 
sample  D-13  is  also  evidence  for  this  model. 

Austin  and  Sheard  ( 1957)  have  measured  the  optical  energy 
gap  of  BigTe^  -  BigSe^  alloys  and  find  a  sharp  break  in  their  curve  at 
BigTe2  ^Be^  If  ve  extrapolate  the  curve  back  to  pure  BigTe^  and 
assume  that  this  energy  difference  of  about  0.2  eV  is  the  energy 
difference  between  the  two  conduction  bands,  we  find  that  the  lower 
hand  has  a  density-of-states  of  0.05mQ.  This  is  Just  the  value  required 
to  explain  the  large  mobility. 


12 


IV.  ACKNOWLEDGMENTS 

We  wish  to  thank  G.  R.  Wagner  for  carefully  performing  most 
of  the  measurements  reported  here,D*A.Zupon  for  assistance  with  crystal 
growing,  R.  C.  Miller  and  A.  Sagar  for  several  stimulating  discussions, 
J.  McHugh  for  growing  some  of  the  crystals,  and  C.  B.  Satterthwaite  for 
the  Seebeck  coefficient  data. 


REFERENCES 

Allgaier,  R.  S.,  and  Scanlon,  W.  W. ,  1958,  Phys.  Rev.  Ill,  1029. 

Austin,  I.  G.,  1958,  Proc.  Phys.  Soc.  ]2,  5*4-5. 

Austin,  I.  G.,  i960,  Proc.  Phys.  Soc.  j6,  169. 

Austin,  I.  G.,  and  Sheard,  A.,  1957,  J-  Electronics  and  Control  21*46. 

Bowley,  A.  E.,  Delves,  R.,  and  Goldsmid,  H.  J.,  1958,  Proc.  Phys.  Soc.  J2,  *401. 
Brooks,  H. ,  1955,  Advances  in  Electronics  and  Electron  Physics,  edited  by 
L.  Marton  (New  York:  Academic  Press). 

Chandrasekhar,  B.  S.,  1961,  Rev.  Sci.  Instr.  £2,  368. 

Cole,  M.,  Bucklow,  I.  A.,  and  Grigson,  C.  W.  B.,  1961,  Brit.  J.  Appl.  Phys. 

12,  296. 

Dingle,  R.  B.,  1955,  Phil.  Mag.  46,  831. 

Drabble,  J.  R. ,  Groves,  R.  D.,  and  Wolfe,  R.,  1958,  Proc.  Phys.  Soc.  J1 ,  *430. 
Goldsmid,  H.  J. ,  1958,  Proc.  Phys.  Soc.  Jl,  633. 

Goldsmid,  H.  J. ,  1961,  J.  Appl.  Phys.  32,  2198. 


13 


REIERENCES  ( continued ) 

Hall,  E.  N.,  i960,  Proceedings  of  the  Intematioial  Conference  on  Semiconductor 
Physics  (Prague:  Czechoslovak  Academy  of  Sciences). 

Herring,  C.,  i960,  Proceedings  of  the  International  Conference  on  Semiconductor 
Physics  (Prague:  Czechoslovak  Academy  of  Science),  p.  63* 

Kanai,  Y.,  Nii,  R.,  and  Watanahe,  N.,  1961,  J.  Appl.  Phys.  j|2,  2146. 

Mansfield,  R.,  and  Williams,  W.,  1958,  Proc.  Phys.  Soc.  ]2,  733* 

Offergeld,  G.,  and  Van  Cakenberghe,  J. ,  1959,  J*  Phys.  Chem.  Solids  11,  310. 

Satterthvreiite,  C.  B.  and  Ure,  R.  W.,  Jr.,  1957*  Phys.  Rev.  108,  1164. 

Sladek,  R.  J. ,  1959,  private  communication;  reported  by  Ure,  1961. 

Ure,  R.  W.,  Jr.,  1961,  Thermoelectricity:  Science  and  Engineering,  edited  by 

R.  R.  Heikes  and  R.  W.  Ure  (interscience  Publishers,  New  York,  196l),p.  413. 

Vinogradova,  M.  N.,  Golikova,  0.  A.,  Efimova,  B.  A.,  Kutasov,  V.  A,, 

Stavitskaya,  T.  S.,  Stil'bans,  L.  S.,  and  Sysoeva,  L.  M. ,  1959* 

Fiz.  Tverdogo  Tela  1,  1333  [English  translation:  Soviet  Phys. -Solid 
State  1,  1224  (i960)]. 

Walker,  P.  A.,  i960,  Proc.  Phys.  Soc.  ]6_,  113* 

Yates,  B.,  1959,  J*  Electronics  and  Control,  6,  26. 


14 


Table  I 

Temperature  Dependence  of  Mobility  in 

N-Type  Tellurium  Doped  Bismuth  Telluride 

Carrier 

Temperature 

Concentration 

Dependence 

Sample 

(cm  3) 

of  Mobility 

D-13 

2.4  x  1017 

x-e 

M3-1 

5-3 

T-2.7 

El-2 

3.0  x  1018 

T-2.2 

A2-1 

3.4 

t-2.4 

A 4-1-5 

1.2  x  1019 

t-1*70 

IU-1 

6.8 

t-1.31 

Table  II 

Values  for  Denaity-of-States  Effective  Mass 

Derived 

from  Seebeck  Coefficient  Data 

at  100°K 

Sample  Number 

D-13  27 

16 

15  10 

Source  of  Data  This  paper  - — — 

. -Walker  (i960)— . 

"Carrier  Concentration"  2.4  x  10^  1.6  x  10l8 

6  x  1018 

1.6  x  1019  4  x  1019 

m,/in  for  s  = 
a'  o 

-1/2  0.26  0.64 

O.67 

0.96  1.00 

m,/m  for  s  = 
a/  o 

3/2  .07  0.21 

0.21 

0.31  0.34 
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Fig.  1 


Fig.  2 


Fig.  3 


Fig.  h 


Fig.  5 


FIGURE  CAPTIONS 

Hall  coefficient  for  tellurium  doped  n-type  bismuth  telluride. 

The  symbols  are  identified  on  Fig.  2.  The  Hall  coefficient  is 
measured  with  the  current  in  the  cleavage  plane  and  the  magnetic 
field  perpendicular  to  the  cleavage  plane. 

Electron  Hall  mobility  in  tellurium  doped  n-type  bismuth  telluride. 
Samples  designated  by  solid  symbols  are  more  homogeneous,  i.e., 
the  two  Hall  coefficients  differed  by  less  than  19j&  and  the  two 
resistivities  differed  by  less  than  10J&. 

Electron  Hall  mobility  vs  "electron  concentration"  at  k.2°K  and 
77°K.  Yates  (1959)  measurements  at  4.2°K  on  iodine  doped  samples 
are  designated  by  a  cross.  All  other  points  are  our  measurements 
on  tellurium  doped  samples.  Samples  designated  by  solid  symbols 
are  more  homogeneous  than  those  designated  by  open  or  partially 
open  symbols. 

Electrical  resistivity  for  tellurium  doped  n-type  bismuth 
telluride.  The  symbols  are  identified  on  Fig.  2. 

Seebeck  coefficient  for  n-type  tellurium  doped  and  p-type  zone 
refined  bismuth  telluride. 
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Fig.  I  -  Hal!  coefficient  for  Te  doped  Bi2Te3 
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Fig,  2  -  Electron  Hall  mobility  for  Te  doped  Bi2Tej 
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Fig.  4  -  Resistivity  of  Te  eloped  Bi4Te5 
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ABSTRACT 

The  absolute  Seebeck  coefficient,  electrical  resistivity, 
and  thermal  resistivity  vere  simultaneously  measured  on  pure  bismuth 
single  crystals  of  various  orientations  between  approximately  80°K  and 
300°K.  Using  an  overlapping  two-band  many- valley  model,  numerical  values 
for  the  temperature  dependence  and  anisotropy  (where  appropriate)  of  the 
following  parameters  have  been  calculated:  (l)  the  overlap  energy  and  the 
Fermi  energy  of  the  electrons  and  of  the  holes,  (2)  the  density  of  states 
effective  mass  of  the  electron  and  of  the  holes,  (3)  the  separate  electronic 
and  lattice  thermal  conductivities,  (4)  the  actual  index  of  thermoelectric 
efficiency,  and  the  hypothetical  "optimum"  index  of  thermoelectric 
efficiency.  The  calculated  electronic  thermal  conductivity  includes 


a  new  term  aue  to  oipolar  diffusion. 


INTRODUCTION 


Bismuth  is  a  semimetal  which  has  a  crystal  structure  based 
on  the  rhombohedral  lattice  with  two  atoms  permit  cell.  Many  of  the 
electrical  and  thermal  properties  of  bismuth  are  anisotropic  and  can 
be  described  with  reference  to  the  trigonal  axis  and  three  equivalent 
binary  axes  lying  in  a  plane  perpendicular  to  the  trigonal  axis. 

The  Jones '  zone  for  bismuth  can accommodate  exactly  five 
electrons  per  atom,  so  that  there  are  just  enough  valence  electrons 
to  fill  a  band.  However,  there  is  an  overlap  in  energy  between  the 
electron  and  hole  bands  and  a  small  fraction  of  the  electrons  spill 
over  into  a  higher  zone,  even  at  0°K.  Information  concerning  the  band 
structure  of  bismuth  has  been  obtained  from  measurements  of  the  de  Haas- 
van  Alphen  effect,1-^  cyclotron  resonance, ^  galvanomagnetic  effects, 10-12 
Seebeck  coefficient,1^  piezoresistance, ^  magneto-optical  effects,1^*1 

117  lO 

ultrasonic  attenuation  and  anomalous  skin  effect. 

We  have  analyzed  our  measurements  in  terns  of  a  parabolic 
two-band  model.  For  both  the  electrons  and  holes,  the  constant  energy 
surfaces  are  ellipsoids  vhich  are  arranged  in  momentum  space  in  accordance 
with  the  rhombohedral  symmetry  of  the  bismuth  crystal.  In  our  calculations, 
we  have  used  Abeles  and  Meiboom's1^  results  for  the  number  of  carriers  and 
the  mobility  ratios.  They  used  the  following  parabolic  two-band  model 
which  is  more  specific  than  ours.  The  energy  surfaces  of  the  holes  are 
spheroidal  with  the  axis  of  revolution  parallel  to  the  trigonal  axis  of 
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the  crystal.  The  energy  surfaces  of  the  electrons  are  ellipsoids,  each 

one  having  one  of  its  axes  parallel  to  the  trigonal  axis  and  another  axis 

l_k 

parallel  to  the  binary  axis.  De  Haas-van  Alphen  effect,  cyclotron 
7— Q  T7  lft 

resonance'”''  and  other  measurements  indicate  that  the  electron 
ellipsoids  are  actually  slightly  tilted  away  from  the  trigonal  plane. 
However,  this  should  have  a  negligible  effect  on  our  results. 
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EXPERIMENTAL 

2 

Single  crystals  of  about  0.2  cm  uniform  cross-section  and 

over  10  cm  long  were  grown  in  a  graphite  mold  in  vacuum,  using  a 

13 

modified  Bridgman  technique  as  previously  described.  The  bismuth, 

of  99*999$  purity,  was  obtained  from  the  American  Smelting  and  Refining 

13 

Company.  Previous  measurements  indicated  that  the  Seebeck  coefficient 
was  essentially  unaffected  by  zone  refining,  thus  indicating  a  high  degree 
of  purity.  Measurements  on  some  crystals,  as  grown,  indicated  a  ratio 
of  room-temperature  to  liquid-helium- temperature  resistivities  of  approxi¬ 
mately  ninety.  Orientations  close  to  the -trigonal  axis  were  obtained  by 

suitable  seeding.  The  orientations  of  the  crystals  were  determined  by 

19 

back-reflection  Laue  photographs  and  further  checked  by  measuring  the 
Seebeck  coefficient  (a)  at  300°K  and  using  a  previously  published^  graph 
of  a  vs.  orientation. 

The  samples  were  carefully  cut  to  approximately  1  cm  length 
with  a  jeweler's  saw  and  the  ends  were  lapped  flat  and  parallel.  The 
samples  were  etched  in  nitric  acid  and  the  ends  were  nickel  plated. 

A  low  melting  point  alloy  (Cerrolov  117)  was  used  to  solder  the  sample 
between  the  sample  heater  and  cold  plug  of  the  apparatus  shown  in  Pig.  1. 

The  apparatus  of  Fig.  1  was  used  to  measure  the  electrical 
resistivity,  thermal  resistivity  t  and  Seebeck  coefficient  of  the 
sample.  The  temperatures  were  measured  with  copper-constantan  thermo¬ 
couples  which  were  thermally  bonded  to  isothermal  regions  before  being 
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brought  out  of  the  apparatus.  All  of  the  measurements  were  made  with 

a  Leeds  and  Northrup  K-3  potentiometer  and  appropriate  standard  resistors. 

The  isothermal  electrical  resistivity  was  measured  by  passing 

a  dc  current  of  about  0.1  amp  through  the  sample.  A  reversing  procedure 

was  used  to  eliminate  thermoelectric  effects.  The  potential  drop  across 

the  sample  was  measured  across  the  copper  leads  of  the  thermocouples  in 

the  sample  heater  and  cold  plug.  Thus,  the  measured  resistance  included 

the  contact  resistance  of  the  soldered  joints.  In  order  to  demonstrate 

that  this  contact  resistance  was  small,  tests  were  made  on  the  assembled 

sample  heater,  sample, and  cold  plug  before  they  were  placed  in  the 

apparatus  of  Fig.  1.  Room  temperature  resistivity  measurements  with 

20 

a  standard  4-probe  ac  apparatus  indicated  that  the  contact  resistance 

and  inhomogeneity  effects  were  less  than  2$.  UBing  an  apparatus  pre- 
21 

viously  described,  potential  profile  measurements  were  also  made  on 
the  assembled  unit  at  room  temperature.  The  linear  nature  of  the  curves 
indicated  that  the  contact  resistance  was  less  than  2$  of  the  sample 
resistance)  that  there  were  no  cracks  in  the  sample)  and  that  the 
electrical  properties  were  uniform  along  the  length  of  the  sample. 

The  Seebeck  coefficient  and  thermal  resistivity  were  measured 
by  establishing  a  steady  state  temperature  gradient  across  the  sample 
and  thermally  matching  the  shield  shown  in  Fig.  1.  The  Seebeck  voltage 
was  measured  across  the  copper  leads  of  the  thermocouples.  The  absolute 
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Seebeck  coefficient  of  bismuth  was  found  by  making  a  correction  for 

22 

the  small  absolute  Seebeck  coefficient  of  copper.  Simultaneously, 

the  thermal  resistivity  was  determined  by  measuring  the  power  input 

21 

to  the  sample  heater.  The  stray  heat  losses  from  the  sample  heater 
(mostly  radiation)  were  determined  as  a  function  of  temperature  in  a 
separate  calibration  run  and  appropriate  corrections  were  made.  The 
total  correction,  which  included  small  corrections  for  the  thermal 
contact  resistance2-^  (<  1$)  and  small  temperature  drifts2^  (<  1.9jt), 
never  exceeded  6$. 

In  order  to  verify  that  the  effects  of  cold  work  were  small 
in  this  temperature  range,  one  sample  was  measured  in  the  usual  way  and 
then  was  severely  strained  by  bending.  Remeasurement  showed  that  the 
properties  had  changed  by  about  &fi  at  100°K. 
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RESULTS  AND  DISCUSSION 

The  absolute  Seebeck  coefficient  (a)>  electrical  resistivity 
(p)  and  thermal  resistivity  (r)  were  simultaneously  measured  between 
approximately  8o°K  and  300°K  on  six  pure  bismuth  single  crystals  of 
various  orientations.  All  of  these  properties  are  anisotropic.  The 
rather  large  anisotropy  of  the  thermal  and  electrical  properties  of 
bismuth  is  somewhat  surprising  since  its  crystal  structure  can  be 
considered  to  be  rhombohedral  with  only  a  slight  distortion- from  cubic. 
Although  the  electrical  and  thermal  resistivities  are  tensor  quantities, 
it  is  convenient  here  to  use  the  term  "resistivity"  to  mean  the  measured 
resistance  times  area/length  and  "Seebeck  coefficient"  to  mean  the 
Seebeck  voltage  divided  by  the  temperature  difference.  Then  for  the 

long  thin  rod  geometry  used  here,  both  the  electrical  and  thermal 

2k 

resistivity  obey  the  Thomuon-Voigt  relation  as  one  would  expect  from 
the  uniaxial  symmetry  associated  with  the  rhombohedral  structure: 

.  Rj^  +  (H  |  |  -  Rj)cos2$  (1) 

where  R  is  either  the  thermal  or  electrical  resistivity  and  ©  is  the 

angle  between  the  sample  axis  and  the  trigonal  axis.  The  subscripts 

| |  and  J.  denote  directions  parallel  and  perpendicular  to  the  trigonal 

axis.  Due  to  the  anisotropy  of  the  thermal  conductivity,  the  absolute 

Seebeck  coefficient  shows  a  deviation  from  the  Thomson-Voigt  relation 

13  25 

and  obeys  the  equation  * 

x(a  1 1  -  aj^os2< 

ou  -  ai  + - 5 — 

v  -*■  1  +  (x-l)cos  © 


(2) 
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where  X  =  Kj/k|  | ,  the  ratio  of  the  thermal  conductivities  perpendicular 

and  parallel  to  the  trigonal  axis  respectively.  The  measured  values  of 

and  r^  have  been  used  to  solve  Eq.  1  by  the  method  of  leaBt  squares  to 

obtain  pj  |,  pj,  r| |  and  r^  as  a  function  of  temperature  as  shown  in  Figs. 

2  and  3>  Likewise,  the  measured  values  of  a ^  have  been  used  to  solve  Eq.2 

to  obtain  ot|  |  and  ct|  as  shown  in  Fig.  4.  The  least  squares  analysis  for 

p,  r  and  a  gave  a  "standard  error  of  estimate"  which  never  exceeded  3$. 

All  measurements  are  believed  to  be  accurate  to  within  +3$.  The  Seebeck 

coefficient  measurements  are  in  excellent  agreement  with  those  of 

Chandrasekhar.1^  The  electrical  resistivity  measurements  are  in  good 

lO  26  27 

agreement  with  those  of  other  authors  '  *  who  used  methods  which 

eliminated  the  electrical  contact  resistance  problem.  The  thermal  con- 

27 

ductivity  measurements  are  in  good  agreement  with  those  of  other  authors, 
although  the  temperature  dependence  of  the  present  results  appears  to  be 
somewhat  slower  at  low  temperatures. 

In  the  following  sections,  several  interesting  parameters  will 
be  calculated: 

(a)  The  Overlap  Energy,  Fermi  Energies  and  Density  of  States  Effective  Mass' 

iq 

It  has  been  previously  shown  that 


ae\ae  * 
Ue}i  + 


(3) 
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where  i  stands  for  | |  or  J_,  a  is  the  electrical  conductivity  and  the 

subscripts  e  and  h  refer  to  electrons  and  holes  respectively.  If  the 

energy  dependence  of  the  relaxation  time  is  independent  of  direction, 

then  the  partial  Seebeck  coefficients,  and  0£^,  are  independent  of 

direction.  The  anisotropy  of  a  arises  from  the  anisotropic  and  unequal 

carrier  mobilities.  Assuming  parabolic  bands,  ellipsoidal  energy  surfaces, 

8  28 

and  that  the  energy  dependence  of  the  relaxation  time  is  E  ,  then 


where  the  Fermi  integrals  are  given  by 

and  the  reduced  Fermi  energies  are 


ie  -  Ef/KT 


(T) 
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where  E^  and  are  the  Fermi  energies  of  the  electrons  and  holes  measured 

from  the  edges  of  choir  respective  hands.  The  overlap  energy,  Eq,  iB 

measured  from  the  top  of  the  valence  band  to  the  bottom  of  the  conduction 

band.  For  intra-valley  acoustical  mode  lattice  scattering,  which  is 

assumed  in  these  calculations,  s  is  -  1/2.  Abeles  and  Meiboom^  have 

calculated  the  anisotropic  mobilities  of  the  two  types  of  carriers  at 

80°K  and  300°K  as  s  no  wit  in  Table  1.  From  their  data,  the  mobility  ratios 

(ue/Hv)  appear  to  be  very  weakly  temperature  dependent.  It  was  assumed 

that  the  individual  mobilities  parallel  and  perpendicular  to  the  trigonal 

axis  had  the  form  p  -  cTJ.  The  results  are  very  insensitive  to  this 

assumption.  Using  Abeles  and  Meiboom's  data,  we  solved  this  equation 
* 

for  the  constants  c  and  y  and  the  mobility  ratios  were  calculated  at 
intermediate  temperatures.  With  these  values  and  the  value  of  ct| |  and 
Op  Eq.  3  was  solved  for  ctg  and  ot^  as  functions  of  temperature  as  shown 
in  Fig.  5.  The  values  of  (XQ  and  0^  were  inserted  into  Eqs.  k  and  5  and 
|e  and  ^  were  obtained  as  a  function  of  temperature  by  the  use  of  tables 
of  Fermi  integrals.  Using  these  values,  we  calculated  the  Fermi  energies 
and  the  overlap  energy  from  Eq.  7.  The  results  are  shown  in  Fig.  6.  The 
equations  and  sign  conventions  used  here  are  the  same  as  those  commonly 
used  in  the  treatment  of  semiconductors.  In  our  notation,  EQ  should  be 

* 


For  both  the  electrons  and  holes,  values  of  y  of  approximately-fi.l 
were  found.  The  significance  of  this  will  be  discussed  later  along 
with  the  scattering  mechanism. 
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a  negative  number  for  the  case  of  an  overlap  of  the  valence  and  conduction 

bands.  The  negative  value  found  here  is  an  independent  determination  that 

29  30 

the  conduction  and  valence  bands  overlap  in  bismuth.  In  other  materials  * 

where  there  is  a  question  as  to  whether  there  is  an  energy  gap  or  energy 

overlap,  this  method  might  be  vised  to  give  a  definite  answer.  The  overlap 

energy  is  a  sensitive  function  of  temperature  as  shown  in  Fig.  6.  This  is 

not  unexpected  since  Eq  depends  on  the  lattice  spacing  which  is  temperature 

sensitive  through  thermal  expansion.  Also,  it  is  well  known  that  the  energy 

gap  in  semiconductors  is  generally  a  function  of  temperature. 

The  total  number  (N)  of  electrons  (or  holes)  per  unit  volume 

has  been  calculated  by  Abeles  and  Meiboom10  at  80°K  and  300°K  as  shown 

28 

in  Table  l.  Assuming  that  the  bands  are  parabolic 

N  =  (^)(2moteT)3/2(m#T)3/2Fl/2(s)  (8) 

* 

where  m  T  is  defined  by  Eq.  8  and  may  be  called  the  "total  density  of 
states  effective  mass"  in  units  of  mQ.  This  parameter,  which  is  isotropic 
and  independent  of  the  number  of  ellipsoids,  has  been  calculated  at  80°K 
and  300°K  (Tables  2  and  3)  using  the  values  of  {  and  N.  This  parameter 
will  be  of  interest  in  discussing  the  thermoelectric  efficiency  later  on. 
Similarly,  the  density  of  states  effective  masses  "per  ellipsoid"  for  the 
electrons  and  holes  were  calculated  at  8o°K  and  300°K  (Tables  2  and  3) 
from  the  equation 
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(|)  =  (^)(2mokT)3/V)3/2Fl/2(6)  (9) 

where  n  Is  the  number  of  ellipsoids,  N/n  is  the  number  of  electrons  or 

* 

holes  per  ellipsoid  per  unit  volume,  and  m  is  the  density  of  states 

effective  mass  per  ellipsoid  in  units  of  mQ.  In  this  approach,  n  is 

* 

a  variable  interger.  This  m  corresponds  to  the  effective  mass  which 
is  derived  from  the  de  Haas-van  Alphen  effect,  cyclotron  resonance  and 
other  methods.  These  measurements  determine  the  components  of  the  reci¬ 
procal  mass  tensor  of  an  individual  ellipsoid,  but  do  not  readily  yield 
information  concerning  the  number  of  ellipsoids. 

Electrons:  The  Fermi  energy  of  the  electrons  has  been  determined 
1-^ 

by  Shoenberg  J from  de  Haas-van  Alphen  effect  measurements.  The  present 
results  (Fig.  6)  are  in  excellent  agreement  with  his  value  of  =  0.0177  ev 
at  h°K.  The  density-of-states  effective  mass  of  the  electrons  per  ellipsoid 
has  been  calculated  from  the  components  of  the  effective  mass  tensor  as 
determined  by  other  metho'd^ 1  of  these  other  measurements  were  made 
at  liquid  helium  temperatures.  Table  2  presents  the  results;  where 
necessary,  the  value  of  E^,  was  assumed  to  be  0.0177  ev.  The  density-of- 
states  effective  mass  per  ellipsoid  that  is  determined  by  the  present  work 
depends  on  whether  a  3  or  6  ellipsoid  model  is  assumed.  As  can  be  seen 
from  Table  2,  both  values  are  in  reasonable  quantitative  agreement  with 
the  other  measurements.  We  are  unable  to  choose  between  3  or  6  ellipsoids 
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because  of  the  divergence  of  the  other  results.  From  the  present  results 

(Table  2),  it  appears  that  the  density-of-state6  effective  mass  of  the 

electrons  is  independent  of  temperature  between  80°K  and  300°K.  This 

31 

indicates  that  the  electron  band  is  parabolic. 

Holes :  The  Fermi  energy  of  the  holes  has  been  determined  by 

Brandt^  from  de  Haas -van  Alphen  effect  measurements.  Our  low  temperature 

results  agree  withhis  value  of  E^  =  O.OI56  ev.  Brandt^  has  also  found 
* 

m  h  =  0.115  which  is  in  excellent  agreement  with  our  low  temperature 

value  that  was  calculated  assuming  one  hole  spheroid  (Table  3).  This 

agreement  is  evidence  that  the  energy  surface  for  the  holes  is  described 

32 

by  only  one  spheroid  .  From  anomalous  skin  effect  measurements  and 

l8 

with  the  assumption  that  there  are  two  hole  spheroids,  Smith  has  found 

values  of  the  components  of  the  reciprocal  mass  tensor  divided  by  the 

Fermi  energy.  Using  Brandt's  value  of  E^  =  O.OI56  ev  with  Smith's  data, 

* 

one  finds  m  .  ~  0.1  m  . 

h  o 

l8  3  3 

From  low  temperature  specific  heat  measurements,  estimates 

of  m^  3f  1  rr.Q  and  5 ,>001  ev  have  been  made.  Recent  galvanomagnetic 

3k 

and  thermomagnetic  experiments  by  Sybert,  Grenier  and  Reynolds  indicate 
the  existence  of  two  hole  bands.  The  uppermost  valence  band  has  a  high 
mobility  and  presumably  a  low  effective  mass.  The  lower  valence  band 
has  a  lower  mobility  and  presumably  a  higher  effective  mass.  Thus  it 
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seems  reasonable  that  the  "dominant"  hole  band  has  a  light  effective 

I  7  1  o  £  *1  Q 

mass  as  found  by  Galt,  et  al.,  Chandrasekhar,  J  Brandt0  and  Smith 

j  (as  modified  above).  Presumably  then,  it  is  this  lower  heavy-mass  band 

which  contributes  to  the  specific  heat.  Also,  this  additional  band  may 

* 

|  be  the  cause  of  the  temperature  dependence  of  m  ^  which  we  have  observed. 

From  the  temperature  dependence  of  E^  (Fig.  6),  it  seems  that  this  extra 
hole  band  plays  a  role  in  the  transport  properties  above  approximately 
200°K. 


Scattering  Mechanism:  Our  results  are  sensitive  to  the  assumption 

of  intra-valley  carrier  scattering  by  acoustical  mode  lattice  waves 

(s  =  -  l/2).  Although  they  are  not  presented,  we  have  also  done  the 

calculations  assuming  that  s  =  +  i/2.  However,  the  Fermi  energies  of 

both  the  electrons  and  holes  as  determined  by  the  de  Haas-van  Alphen 
1-’  6 

measurements  are  in  agreement  with  our  results  only  if  a  value  of 

s  —  -  1/2  is  assumed.  Thus,  at  low  temperatures,  it  appears  that  both 
the  electrons  and  holes  are  predominantly  scattered  by  acoustical  mode 

lattice  waves  as  one  would  expect.  For  such  a  case,  theory  predicts 

“1**)  _2«1 
u  -  cT  .  As  previously  described,  however,  we  have  found  p.  S’  cT 

from  Abeles  and  Meiboom's  data.  This  discrepancy  is  probably  due  to  the 
effect,  of  the  second  valence  band.  However,  it  may  also  be  due  to  non¬ 
parabolic  hole  band(s),  the  onset  of  optical  mode  lattice  scattering  or 
inter- valley  effects.  In  any  event,  our  treatment  does  not  take  these 
higher  temperature  effects  into  account. 
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(b)  The  Electronic  and  Lattice  Thermal  Conductivities 

It  was  assumed  that  the  measured  thermal  conductivity  is 
the  sum  of  the  lattice  thermal  conductivity  (K^)  and  the  total  electronic 
thermal  conductivity  (Kg). 

K  =  +  Kg  (10) 

* 

For  the  present  model,  the  total  electronic  thermal  conductivity  is 
given  by  the  sum  of  three  terms:  J  (l)  the  ordinary  thermal  conductivity 
due  to  electrons  alone,  *e,  (2)  the  ordinary  thermal  conductivity  due 
to  holes  alone,  and  (3)  the  transport  of  heat  due  to  bipolar  diffusion,*^. 


*E  =  *e  +  Kh  +  *eh 


(11) 


The  ordinary  thermal  conductivities  have  been  extensively  treated  in 
28 

the  literature.  However,  the  bipolar  diffusion  term  for  the  case  of 

35 

a  semi -metal  with  overlapping  bands  has  only  recently  been  realized  — 
although  the  analogous  situation  in  semiconductors  is  well  known. 
Briefly  this  term  is  present  for  the  following  reason.  The  thermal 
conductivity  is  defined  as  the  heat  flux  divided  by  the  temperature 


As  previously  discussed,  there  really  appear  to  be  more  than  two 
active  bandB  in  bismuth.  Therefore,  there  will  be  other  terms  in 
K_.  However,  these  extra  terms  will  be  small  and  we  have  done  the 
calculations  on  the  two-band  model. 
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gradient  when  the  total  electrical  current  is  zero.  The  ordinary  terms 
give  the  heat  flow  contributions  of  the  electrons  and  holes  separately 
when  their  individual  electrical  currents  are  zero.  For  materials  with 
carriers  of  both  signs,  there  is  an  additional  term  which  arises  because 
the  individual  electrical  currents  are  not  zero  although  the  total  electri¬ 
cal  current  is  zero.  This  effect  may  be  thought  of  as  the  transport  of 
the  "heat  of  formation"  of  electron-hole  pairB  which  diffuse  from  the  hot 
region  to  the  cold.  From  our  previous  calculations,  it  is  possible  to 
evaluate  all  three  terms,  which  are 


7(Se)°e  +  ^h>a 


h 


(12) 


where 


7(6)  = 


(s  4  7/2)Fb  +  ,/2(|) 

■(s  y2[t) 


(B  4  5/2)Fb  +  3/2(t) 

4  73311 


2 


(13) 


The  results  for  the  principal  directions  are  shown  in  Figs.  J,  7  and  8. 

k  2 

The  values  of  (-)  7(5)  that  were  found  from  Eq.  13  and  used 
in  calculating  *e  and  are  shorn  in  Table  4.  These  values  of  Ke/a^S 
and  Kjj/ojjT  are  considerably  different  from  the  "usual  Lorenz  Humber" 
of  2.45  x  10*®(volts/°C)^,  which  only  holds  in  the  degenerate  approxi¬ 
mation.  The  carriers  in  bismuth  must  be  described  by  the  exact  Fermi- 
Dirac  statistics  and,  in  such  a  case,  7(5)  is  a  function  of  the  Feral 
level  and  the  scattering  mechanism  as  is  given  in  Eq.  13.  The  energy 
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dependence  of  the  relaxation  time  has  been  assumed  to  be  Independent  of 
direction,  which  implies  that  7(1)  is  isotropic.  In  Table  4,  it  can  be 
seen  that  r(te)  -  y(lh)'  Thus,  letting  7  -  7({e)  -  7(th)>  Eq.  12  can 
rewritten  as 


(1*0 


In  order  to  compare  our  results  with  earlier  measurements  on  bismuth 
crystals,  the  "pseudo-Lorenz  Humber,"  L,  has  been  calculated  from  our 
data  as  shown  In  Table  4.  ThiB  quantity,  which  presumably  1b  the  one 
determined  experimentally  in  the  work  discussed  below,  is  defined  as 


L 


(15) 


How,  if  the  bipolar  term  were  not  present  in  Eq.  14,  L  would  be 

isotropic.  However,  from  magnetoresistance  and  magneto-thermal 

36 

conductivity  measurements,  Gruneisen  and  Gielessen  have  found  that 
L  is  anisotropic.  The  bipolar  term  (*eh)  provides  such  an  anisotropic 

O 

contribution  to  L  because  the  factor  is  anisotropic.  The 

qualitative  agreement  between  the  present  calculations  and  Gruneisen 
and  Gielessen's  measurements  iB  Btrong  evidence  for  the  existence  of 
the  bipolar  contribution  to  the  electronic  thermal  conductivity. 


IT 

The  magnetoresistance  and  magneto-thermal  conductivity  measure- 
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ments  of  Gruneisen,  Rausch  and  Weiss^  at  60°K  perpendicular  to  the  tri¬ 
gonal  axis,  yield  a  value  of  L  in  reasonable  agreement  with  our results. 

It  is  very  interesting  that  the  value  of  L  found  by  Gruneisen,  Rausch  and 
Weiss,  as  veil  as  some  of  those  calculated  here  (Table  ^),  are  close  to 
the  "visual  Lorenz  Number"  —  a  pure  accident.'  The  agreement  is  merely 
due  to  the  fortuituous  combination  of  the  bipolar  term  and  the  ordinary 

terms  properly  derived  with  exact  Fermi -Dirac  statistics.  The  values  of 

58 

L  obtained  by  de  Haas,  Gerriteen  and  Capel  by  saturation  magneto-thermal 
conductivity  measurements  at  SonK  parallel  to  the  trigonal  axis  are  not  in 
agreement  vith  our  calculations.  However,  their  particular  type  of  measure¬ 
ment  may  not  be  considered  a  satisfactory  method  of  separating  out  tine 
electronic  thermal  conductivity.  In  view  of  the  present  situation,  careful 
experimental  and  theoretical  treatments  of  magnetoresistarice  and  magneto- 
thermal  conductivity  would  be  valuable,  particularly  at  high  fields. 

The  lattice  thermal  conductivity  was  calculated  from  Eq.  10 
for  each  of  the  principal  directions.  The  lattice  thermal  resistivity, 
which  is  the  scalar  reciprocal  of  the  conductivity  for  the  principal 
directions,  is  shown  in  Fig.  9  as  a  function  of  temperature. 

Theoretically,  the  lattice  thermal  conductivity  is  given 
roughly  by 


a  1/3  cvX 


(16) 
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where  c  is  the  lattice  specific  heat  per  unit  volume,  v  is  the  velocity 
of  sound  and  X  ia  the  phonon  mean  free  path.  The  lattice  specific  heat 
of  bismuth  is  approximately  given  by  the  Debye  theory  with  a  Debye  tempera- 
ture  (9)  of  120  K.  The  velocity  of  sound^  is  essentially  independent  of 
temperature  and  it  is  anisotropic,  which  at  least  in  part  accounts  for  the 
anisotropy  of  the  lattice  thermal  conductivity.  For  an  isotropic  solid 
above  the  Debye  temperature,  Ioffe  has  recently  proposed  that  for 
phonon-phonon  Umklapp  scattering  that 

T 

i  =  A^cdT  (IT) 

where  A  is  a  constant  for  a  given  material.  The  phonon  mean  free  path 
at  each  temperature  was  roughly  calculated  from  Eq.  16  for  the  principal 
directions.  The  temperature  dependence  of  X  obtained  in  this  way  was 
found  to  follow  Eq.  17  within  the  rather  poor  precision  of  the  comparison. 
Above  the  Debye  9,  the  temperature  dependence  of  the  lattice  thermal 
resistivity  is  essentially  contained  in  the  l/x  factor.  At  these  tempera¬ 
tures,  c  is  nearly  constant  and  it  can  be  shown  that  the  integral  in 
Eq.  17  is  proportional  to  (T  -  9/5).  Thus,  above  the  Debye  9,  Ioffe**0 
proposes  that  the  lattice  thermal  resistivity  should  approach 

rL  CC(T  -  9/5)  . 


(18) 
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According  to  Eq.  18,  a  plot  of  lattice  thermal  resistivity  versus 
absolute  temperature  should  approach  a  straight  line  which  should 
extrapolate  back  to  T  =  6/3 •  From  Fig.  9  ,  it  can  be  seen  that  our 
results  are  in  better  agreement,  with  Eq.  18  than  with  older  theories 
which  predict  that  for  T  >  6,  the  lattice  thermal  resistivity  should 
extrapolate  to  zero  at  0°K.  However,  the  agreement  with  Ioffe's  theory 
is  only  qualitative.  As  one  would  expect  from  Ioffe's  Ideas,  there  is 
a  deviation  from  Eq.  18  in  Fig.  9  around  the  Debye  temperature. 


(c) 


Thermoelectric  Efficiency 

28 

The  thermoelectric  "index  of  efficiency"  (ZT)  is  given  by 


ZT  =  (^)T 


(19) 


This  parameter,  which  has  been  calculated  for  the  two  principal 

orientations  from  our  data,  is  shown  in  Fig.  10.  For  real  bismuth, 

1*1 

the  maximum  value  of  Z  occurs  parallel  to  the  trigonal  axis. 

Since  the  macroscopic  parameter  Z  is  a  sensitive  function 
of  the  carrier  concentrations,  it  is  convenient  to  speak  of  a 
"hypothetical  optimum  index  of  efficiency"  (ZqT).  Briefly,  ZQ  occurs 
when  only  one  type  of  carrier  is  present  and  the  Fermi  level  is 
adjusted  to  yield  an  optimum  value  of  Z.  The  presence  of  only  one 
type  of  carrier  ensures  that  the  Seebeck  coefficient  (Eq.  3)  is  a 
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maximum  and  that  thermal  conduction  by  bipolar  diffusion  is  essentially 

26 

zero.  It  can  be  shown  that  ZqT  is  a  monotonically  increasing  function 
of  the  microscopic  material  parameter  (p)  defined  as 


where  nc  is  the  classical  carrier  mobility,  which  is  related  to  the 
actual  mobility  by  the  relation 

uc.  =  n 

From  the  data  previously  discussed  and  Eqs.  20  and  21,  p  has  been  calculated 
assuming  that  nc,  m  T,  and  are  unchanged  in  meeting  the  physical  require¬ 
ments  mentioned  above.  Knowing  P,  the  hypothetical  "optimum"  index  of 

pQ 

efficiency  (ZqT)  was  obtained  from  graphs  published  by  Ure.  The  results 
for  the  two  principal  orientations  are  shown  in  Fig.  11  for  N-type  material. 
The  results  for  P-type  material  do  not  look  promising  and  are  not  presented. 

The  anisotropy  of  ZqT  (~- 1.6:1)  arises  from  the  anisotropic 
nature  of  the  ratio  Ug/K^.  In  addition  to  the  anisotropy  of  UgAj^  the ' 
larger  anisotropy  of  ZT  (—  6.5:1)  in  real  bismuth  is  due  to  the  anisotropy 
of  the  absolute  Seebeck  coefficient  and  thermal  conduction  by  bipolar 
diffusion.  The  anisotropy  in  a  and  is,  in  turn,  due  to  the  aniso¬ 
tropy  in  the  mobility  ratios,  u  /p.. 

v  U 


(s  +  1/2).'F1/2(S) 
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The  Z  of  real  bismuth  Is  considerably  below  ZQ  primarily  because 

of  the  presence  of  the  holes  which  adversely  effect  a  and  Kg.  Also,  the 

electron  concentration  is  too  high  to  give  the  optimum  Z.  By  suitable 

alloying  with  antimony  (e»10)6),  the  bands  can  be  separated  and  a  discrete 

k2 

energy  gap  introduced.  With  appropriate  N-type  doping,  the  deleterious 

effect  of  the  holes  can  be  minimized  and  the  electron  concentration  can  be 

adjusted  to  yield  an  optimum  Z.  Experiments  of  this  type  by  Smith  and 
k^ 

Wolfe  J  indicate  that  reasonably  high  values  of  ZT  can  be  obtained.  In 

fact,  it  is  conceivable  that  a  Bi-Sb  alloy  with  suitable  doping  may  yield 

a  value  of  ZT  higher  than  the  ZqT  shown  in  Fig.  11.  Due  to  the  difference 

in  the  phonon  and  electron  wavelengths,  scattering  by  Sb  atoms  may  lower 

,  kk 

more  than  ue  and  thus  increase  the  ratio  and  ZqT.  This  argument 

will  be  particularly  true  at  low  temperatures,  and  indeed,  this  appears 

to  be  the  case  experimentally. 

k5  .  • 

It  has  been  argued  that  anisotropic  semiconductors  may  be 

apromislng  area  for  the  search  of  thermoelectric  materials  in  spite  of 
the  fact  that  the  high  mobility  direction  will  probably  coincide  with 
the  high  lattice  thermal  conductivity  direction.  In  bismuth,  at  least, 
the  high  electron  mobility  direction  coincides  with  the  low  lattice 
thermal  conductivity  direction.  This  fact  makes  the  search  for  good 
thermoelectric  materials  among  anisotropic  materials  seem  all  the 
more  promising. 
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SUMMARY 

Our  measurements  have  been  analyzed  In  terms  of  a  parabolic 
two-band  multi-valley  model,  where  the  two  bands  overlap.  However, 
other  recent  work12*^  indicates  that  there  really  are  two  hole  bands. 
This,  together  with  the  present  results,  indicates  that: 

(1)  The  electrons  are  described  by  one  band  with  either 
three  or  six  ellipsoids.  The  electron  band  is  parabolic. 

(2)  There  are  two  non-degenerate  hole  bands.  The  higher 
band,  which  is  dominant  at  low  temperatures,  is  described  by  a' single 
spheroid  with  a  high  mobility,10  a  low  density  of  states  effective 
mass^,^,'L^  (m^  £0.1  mQ),  and  a  Fermi  energy^  comparable  to  that  of 
the  electrons (E^  =0.0l6  ev  at  80°K).  The  lower  hole  band^  has  a 
lower  mobility,  a  higher  density  of  states  effective  masB,  and  a 
lower  Fermi  energy. 

(3)  The  magnitude  of  the  energy  overlap  between  the  electron 
band  and  the  highest  hole  band  is  an  increasing  function  of  temperature. 
Thus  the  Fermi  energies  are  a  function  of  temperature. 

(H)  The  electrons  and  holes  are  predominantly  scattered  by 
acoustical  mode  lattice  waves  at  low  temperatures. 

(5)  The  transport  properties  of  bismuth  can  be  treated  with 
the  equations  developed  for  semiconductors,  provided  that  exact  Fermi- 
Dirac  statistics  are  used.  There  are  undoubtedly  other  semi -metals 
which  can  be  successfully  treated  with  the  present  methods. 
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(6)  In  bismuth,  there  is  a  sizeable  additional  term  in  the 
electronic  thermal  conductivity  due  to  bipolar  diffusion.  Up  till  now, 
this  effect  has  been  overlooked.  Between  100°K  and  300°K,  thermal  con¬ 
duction  by  electrons  and  holes  is  the  dominant  mechanism. 

(7)  The  temperature  dependence  of  the  lattice  thermal 
resistivity  is  better  described  by  the  recent  ideas  of  Ioffe  than 
by  older  theories  which  predict  an  extrapolation  to  zero  at  0°K. 

(8)  From  the  experiments  of  Smith  and  Wolfe^  and  from 
this  study  of  pure  bismuth,  it  appears  that  Bi-Sb  alloys  with  appro¬ 
priate  doping  will,  be  of  practical  thermoelectric  interest.  In  bismuth, 
the  high  electron  mobility  direction  coincides  with  the  low  lattice 
thermal  conductivity  direction.  This  fact  makes  the  search  for  good 
thermoelectric  materials  among  anisotropic  semiconductors  and  semi- 
metals  seem  all  the  more  promising. 
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Table  1.  The  numberandL  mobility  of  electrons  and  holes  in  bismuth 
according  to  Abeles  and  Meiboom.10  N  =  P  is  the  number  of  electrons 

-3  2-1-1 

or  holes  per  cm  .  The  ^ 's  danote  mobilities  in  cm  statvolt  sec  . 
The  subscripts  e  and  h  denote  electrons  and  holes  respectively,  and 
the  subscripts  1 1  and  _[  indicate  directions  parallel  and  perpendicular 
to  the  trigonal  axis  respectively. 


300°K 

8o°K 

N,  P  x  10-18 

2.2 

.U6 

(ue)||  x  lo-6 

5-7 

100. 

(ue)j  x  10'6 

4.87 

85.6 

(uh) | j  x  10’6 

♦  62 

10. 

(uh)j  x  10'6 

2.3 

37- 

MV  1  1 

9.19 

10.0 

2.12 


2.31 
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Table  4.  The  ratio  of  the  carrier  contributions  to  the  thermal  aril 

electrical  conductivities.  Equation  13  vas  used  to  calculate 
2  2 

Ke/oeT  =  (|)  r(te)  and  ich/ahT  a  (|)  y({h).  Values  for  the  "pseudo- 
Lorenz  number,"  L  =  K^/cT,  were  calculated  from  Eqs.  12  and  13  and 
include  the  bipolar  term  of  the  total  electronic  thermal  conductivity. 


T^) 

CM 

o' 

«)  > 

Id  w 

*h/°hT 

(v/°C)2 

lir$Dn 

(v/°C)2 . 

li '  % 

(v/°C)2  1 

100 

1.85  x  lo"8 

1.84  x  10-8 

2.30  x  LO-8 

2.98  x  10”8 

150 

1-79 

1*75 

2.42 

3.35 

200 

1.78 

1.74 

2.45 

3.44 

250 

1.79 

1.79 

2.38 

3.26 

300 

1.81 

1.88 

2.29 

3.00 
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FIGURE  CAPTICBS 


Fig.  1.  Schematic  diagram  of  the  apparatus  for  measuring 

a,  p  and  K.  The  sample  vas  soldered  between  the  sample  heater  and  the 

cold  plug.  For  Seebeck  coefficient  and  thermal  resistivity  measurements, 

the  temperature  of  the  thermal  shield  vas  matched  to  that  of  the  sample 

-5 

heater.  The  parts  shown  were  suspended  and  evacuated  to  10  mm  Hg  in  a 
can  which  was  placed  in  an  appropriate  low  temperature  hath.  The  overall 
temperature  was  adjusted  with  a  heater  in  the  heat  sink. 

Fig.  2.  Electrical  resistivity  vs.  temperature  for  pure 
bismuth  single  crystals ;  pj  |  and  pj  were  calculated  by  fitting  Eq.  1 
to  the  data. 


Fig.  3*  The  total  measured  thermal  conductivity  (K)  and 
the  total  calculated  electronic  thermal  conductivity  (Kg)  vs.  temperature 
for  pure  biBmuth  single  crystals.  Ki 


by  fitting  Eq.  1  to  the  data. 
Eqs.  11-13. 


(Kg)  |  j  and  JlCg)^ 


-i-  and  Ki  "  —  were  calculated 
r  1 1  X  r  i 


were  calculated  from 


Fig.  4.  Absolute  Seebeck  coefficient  ve.  temperature  for 
pure  bismuth  single  crystals;  a  1 1  and  a|  were  calculated  by  fitting 
Eq.  2  to  the  data. 


33 


FIGURE  CAPTIONS  (continued) 

Fig.  5.  The  isotropic  partial  Seebeck  coefficient  of  the 
of  the 

electrons  (ag)  and^ holes  (a^)  vs.  temperature  for  pure  bismuth  single 
crystals.  The  dotted  lines  are  extrapolated  to  zero  at  0°K. 

Fig.  6.  The  overlap  energy  (Eq)  and  the  Fermi  energy  of  the 
electrons  (E^)  and  holes  (E“)  vs.  temperature  for  pure  bismuth  single 
crystals. 


Fig.  7.  The  individual  contributions  to  the  total  electronic 
thermal  conductivity  (Kg)  for  pure  bismuth  single  crystals  perpendicular 
to  the  trigonal  axis.  The  ordinary  thermal  conductivity  due  to  electrons 
alone  (*e),  the  ordinary  thermal  conductivity  due  to  holes  alone  (k^), 
and  the  thermal  conductivity  due  to  bipolar  diffusion  (i»eh)  were  calculated 
from  Eqs.  11-13.  For  this  orientation,  (ue/uh)j  =  2. 

Fig.  8.  The  individual  contributions  to  the  total  electronic 
thermal  conductivity  (Kg)  for  pure  bismuth  single  crystals  parallel  to 
the  trigonal  axis.  The  ordinary  thermal  conductivity  due  to  electrons 
alone  (*e),  the  ordinary  thermal  conductivity  due  to  holes  alone  (k^) , 
and  the  thermal  conductivity  due  to  bipolar  diffusion  (tcell)  were  cal¬ 
culated  from  EqB.  U-13.  For  this  orientation,  (he/Hjj)  1  |  ^  10. 


FIGURE  CAPTIOUS  (continued) 
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Fig.  9.  The  calculated  lattice  thermal  reaiativity  vb. 
temperature  for  pure  bismuth  single  crystals.  The  scatter  at  the 
higher  temperatures  is  due  to  the  fact  that  the  lattice  thermal 
conductivity  is  the  difference  between  two  large  numbers  in  Eq.  10 
as  can  be  seen  in  Fig.  3. 

Fig.  10.  The  index  of  thermoelectric  efficiency  vb. 
temperature  for  pure  bismuth  single  crystals. 

Fig.  11.  The  hypothetical  optimum  index  of  thermoelectric 
efficiency  vs.  temperature  for  "optimized  N-type  bismuth  crystals." 
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Fig.  1 

Schematic  diagram  of  the  apparatus  for  measuring  a,  p  and  K.  The  sample 
was  soldered  between  the  sample  heater  and  the  cold  plug.  For  Seebeck  coefficient  and 
thermal  resistivity  measurements,  the  temperature  of  the  thermal  shield  was  matched  to 
that  of  the  sample  heater.  The  parts  shown  were  suspended  and  evacuated  to  IxT'W  Hg 
in  a  can  which  was  placed  in  an  appropriate  low  temperature  bath.  The  overall  temperature 
was  adjusted  with  a  heater  in  the  heat  sink. 
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Fig.  2 


Electrical  resistivity  vs.  temperature  for  pure  bismuth 
single  crystals;  pit  and.  p>  were  calculated  by  fitting  Eq.  1  to  the  data. 


Thermal  Conductivity  (w/cm-°C) 


Temperature  (  K) 


The  total  measured,  thermal  conductivity  (K)  and  the  total 
calculated  electronic  thermal  conductivity  (Kg,)  vs.  temperature  for  pure 
bismuth  single  crystals.  K|  |  =  an<l  Kj  =  ^  were  calculated  by  fitting 
Eq.  1  to  the  data.  (K_)  i  i  and  o%  ,  were  calculated  from  Eqs.  11-13. 


The  Isotropic  partial  Seeheck  coefficient  of  the  electrons  (afi)  and  of  the  holes  (o^)  vs.  temperature 
for  pure  bismuth  single  crystals.  The  dotted  lines  are  extrapolated  to  zero  at  0°K. 
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rhe  calculated  lattice  thernel  resistivity  vs.  temperature  for  pure  bismuth  single  crystals.  The  scatter 
at  the  higher  temperatures  is  due  to  the  fact  that  the  lattice  thermal  conductivity  is  the  difference 
+.«n  lsrire  numbers  in  Ea#  10, as  can  be  seen  in  Fig.  3* 
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The  hypothetical  optim*  index  of  thermoelectric  efficiency  vs.  temperature  for  “optimized  N-type  bismuth  crystal*" 


